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ABSTRACT 



Advection-dominated accretion flows (ADAFs) have a positive Bernoulli parame- 
ter, and are therefore gravitationally unbound. The Newtonian ADAF model has been 
generalized recently to obtain the ADIOS model that includes outflows of energy and 
angular momentum, thereby allowing accretion to proceed self-consistently. However, 
the utilization of a Newtonian gravitational potential limits the ability of this model to 
describe the inner region of the disk, where any relativistic outflows are likely to origi- 
nate. In this paper we modify the ADIOS scenario to incorporate a pseudo-Newtonian 
potential, which approximates the effects of general relativity. The analysis yields a 
unique, self-similar solution for the structure of the coupled disk/wind system. Inter- 
esting features of the new solution include the relativistic character of the outflow in 
the vicinity of the radius of marginal stability, which represents the inner edge of the 
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quasi-Keplerian disk in our model. Hence our self-similar solution may help to explain 
the origin of relativistic jets in active galaxies. At large distances the radial depen- 
dence of the accretion rate approachs the unique form M oc r 1//2 , with an associated 
density variation given by p oc r _1 . This density variation agrees with that implied by 
the dependence of the X-ray hard time lags on the Fourier frequency for a number of 
accreting galactic black hole candidates. While intriguing, the predictions made using 
our self-similar solution need to be confirmed in the future using a detailed model that 
includes a physical description of the energization mechanism that drives the outflow, 
which is likely to be powered by the shear of the underlying accretion disk. 

Subject headings: accretion — accretion disks — black hole physics — galaxies: jets 

1. INTRODUCTION 

Advection-dominated accretion flows (ADAFs) are attractive models for explaining the emis- 
sion from underluminous black hole candidates such as the Galactic center source. The earliest 
paper on this subject was by Ichimaru (1977), and interest in the subject was revived by Narayan 
& Yi (1994, 1995a, 1995b) and by Abramowicz et al. (1995). Comprehensive reviews of work in 
this area can be found in Kato et al. (1998) and in Narayan, Mahadevan, & Quataert (1998). 
Briefly, these models describe the dynamics of gas fed onto a black hole at very low (significantly 
sub-Eddington) accretion rates. The gas is optically thin and hot, and possesses separate ion 
(Tj ~ 10 12 K) and electron (T e ~ 10 9 K) temperatures because it is too tenuous for Coulomb col- 
lisions to equilibrate the species. The low density also reduces the efficiency of bremsstrahlung 
emission, hence limiting the luminosity of the flows. Although some success has been achieved in 
understanding underluminous systems using ADAF models, a number of issues have been raised 
concerning their self-consistency. Attention has been focused for example on the two-temperature 
assumption, which is one of the basic features of ADAF models. While some authors have identi- 
fied conditions under which this assumption is justified (Blackman 1999; Gruzinov 1998; Quataert 
1998; Quataert & Gruzinov 1999), the validity of the assumption has also been challenged based on 
the possibility that dissipation via plasma instabilities may preferentially heat the electrons rather 
than the protons (Bisnovatyi-Kogan & Lovelace 1997). Chakrabarti (1999) has questioned the rel- 
evance of the self-similar solutions that do not incorporate boundary conditions and therefore do 
not exhibit shocks. 

Narayan, Kato, & Honma (1997) and Narayan & Yi (1994, 1995a) have pointed out that in 
ADAF models the gas has a positive Bernoulli parameter and is therefore gravitationally unbound. 
The positivity of the Bernoulli parameter reflects the fact that the viscous torque transports energy 
as well as angular momentum outward through the disk, which tends to make the total local 
energy per unit mass positive. This could result in the evaporation and escape of at least some of 
the gas before it passes through the event horizon of the black hole. Hence the standard ADAF 
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models may not be self-consistent, since they assume that all of the available gas accretes onto the 
compact object. This idea has been further supported by the ADAF simulations of Igumenshchev 
& Abramowicz (1999), which indicate that bipolar outflows may be produced when the Shakura- 
Sunyaev (1973) viscosity parameter a ~ 1. However, Abramowicz et al. (2000) have suggested that 
self-consistent advection-dominated models with a positive Bernoulli parameter and no outflow can 
be constructed if the self-similar approach is abandoned. Ogilvie (1999) makes a similar point 
and solves a simplified version of one-dimensional, quasi-spherical accretion to conclude that time- 
dependent accretion can take place without the need for any outflows. While this is an interesting 
possibility, it appears that many underluminous black hole systems such as Sgr A* may possess 
relativistic outflows (Falcke 1999; Donea et al. 1999 and references therein). We therefore focus on 
coupled disk/wind models in this paper. 

To address the problem of gravitationally unbound accretion in the standard ADAF models, 
Blandford & Begelman (1999; hereafter BB) generalized the model to include the possibility of 
powerful winds that carry away mass, energy, and angular momentum in proportions that render 
the Bernoulli constant of the gas remaining in the disk negative, thereby allowing it to accrete. 
The resulting advection-dominated inflow-outflow solutions (ADIOS) are similar to the ejection- 
dominated accretion flows (EDAFs) analyzed by Donea et al. (1999). Beckert (2000) has also 
obtained related self-similar solutions for advection-dominated accretion flows that include outflows. 
As it stands, the formalism of BB does not address the critical question of the formation mechanism 
for the outflow. This mechanism must be capable of accelerating high-energy protons, thereby 
creating a nonthermal tail in the particle distribution. Protons in the nonthermal tail escape 
to form the wind whereas the low-energy protons are tightly bound and therefore accrete. As 
discussed by Subramanian, Becker, & Kazanas (1999; hereafter SBK), the energy transfer may be 
accomplished via second-order Fermi acceleration driven by the shearing magnetic field. 

One of the primary reasons for studying black hole accretion is to understand the connection 
between accretion and the production of relativistic jets. The jets probably originate close to the 
black hole as relativistic winds that become collimated with increasing distance. The study of the 
formation of the wind/jet therefore requires a reasonably precise treatment of the inner region of 
the disk, where the physical conditions are most extreme. Although it provides a useful means for 
analyzing steady accretion at large distances from the black hole, the ADIOS model in its original 
form is not applicable close to the event horizon due to the utilization of the standard Newtonian 
form for the gravitational potential. In this paper we attempt to incorporate the effects of general 
relativity by utilizing the pseudo-Newtonian approximation for the gravitational potential suggested 
by Paczyhski & Wiita (1980). Clearly, this approach is not a substitute for a full treatment of 
general relativity, but it does provide us with a convenient mathematical tool for exploring the 
structure of the inner region of the disk. Our strategy is to construct new models for self-similar, 
advection-dominated accretion disks and associated winds that are as faithful as possible to the 
original ADAF and ADIOS scaling principles, while incorporating the pseudo-Newtonian form for 
the gravitational potential. The solutions obtained describe the structure of the entire disk outside 
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the radius of marginal stability, which is the edge of the quasi-Keplerian disk in our model. One 
clear implication of our approach is that our model should reduce to the ADIOS form at large 
distances from the black hole. We focus here on gasdynamical winds that exert no torque on the 
underlying disk. These winds may be powered by the shear acceleration mechanism discussed by 
SBK. The remainder of the paper is organized as follows. In § 2 we develop the pseudo-Newtonian 
analog to the self-similar ADIOS model of BB. The consequences for the structure of the disk/wind 
system are analyzed in §§ 3 and 4, and the implications for the interpretation of observations 
of under luminous black-hole systems are discussed in § 5. This includes a consideration of the 
relevance of the SBK shear acceleration in the context of the pseudo-Newtonian disk/wind model 
developed here. 

2. PSEUDO-NEWTONIAN INFLOW-OUTFLOW MODEL 

Working in units such that GM = c = 1, we can modify the thin-disk, quasi-Keplerian 
ADIOS formalism of BB to approximate the effects of general relativity close to a Schwarzschild 
black hole by replacing the Newtonian potential $ = — l/r with the pseudo-Newtonian form <E> = 
— l/(r — 2) introduced by Paczyhski and Wiita (1980) and further discussed by Paczyhski (1998). 
This prescription yields the correct values for the marginally bound and marginally stable orbital 
radii. The pseudo-Newtonian potential offers a convenient means of approximating the effects of 
general relativity while preserving the degree of mathematical simplicity contained in the ADIOS 
solutions. However, as we demonstrate below, the pseudo-Newtonian potential introduces a length 
scale (namely, the horizon at r = 2), and therefore it is no longer possible to construct self-similar 
models based solely on power-law variations in radius. 

2.1. Self-Similar Forms 

Our goal here is to use the same physical principles underlying the ADIOS scalings to obtain 
new self-similar forms that are valid in the inner region of a disk governed by the pseudo-Newtonian 
potential. The advantage of the self-similar approach is that the solutions obtained may be valid 
over several orders of magnitude in radius. As with any self-similar model, there is a certain degree 
of arbitrariness in the imposed scalings. We therefore base our approach on the same physical 
principles employed in the ADIOS model, as applied in the context of the pseudo-Newtonian po- 
tential. We believe that the self-similar solutions obtained can yield useful insights into possible 
behaviors that should be investigated subsequently using models based on rigorous descriptions of 
the relevant physics. 

BB assumed that the angular velocity of the disk 0, satisfies the quasi-Keplerian scaling f2 oc 
VLk-, where fix oc r~ 3 / 2 is the Keplerian angular velocity in a Newtonian potential. In the pseudo- 
Newtonian case the appropriate expression for the Keplerian angular velocity is (Paczyhski &: Wiita 
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1980) 

-1/2 

and therefore the variation of Q in a quasi-Keplerian accretion disk now becomes 

n(r) = ^ = n r -^, (2.2) 
r r — 2 

where is a positive constant and v$ is the azimuthal velocity. Hence f2 scales as Q.k throughout 
the disk. As in the ADIOS model, we assume that < Uq < 1 so that the disk is sub-Keplerian. 
Simulations carried out by Narayan et al. (1997) indicate that the behavior of 12 is well represented 
by equation (2.2) over several orders of magnitude in radius (outside the radius of marginal stability) 
for values of the Shakura-Sunyaev viscosity parameter a ~ 0.1 — 0.3. Based on these results, we shall 
assume in our work that the variation of 12 is given by equation (2.2). However, this assumption 
does not apply in the region below the radius of marginal stability, where the gas passes through 
a sonic point before crossing the horizon. We provide further discussion of this point in § 5. 

Proceeding by analogy with the ADAF/ADIOS models, we shall suppose that the outwardly- 
directed energy transport rate E and the inwardly-directed angular momentum transport rate / in 
the disk scale as 

E(r) = e^, f(r) = A*, (2.3) 

respectively, where e and A are positive constants, £ = r 2 12 is the angular momentum per unit 
mass, and M(r) is the positive mass accretion rate. We will also assume that the internal energy 
per unit mass (i.e., the ion temperature) scales as the gravitational potential energy, implying that 

where a is the isothermal sound speed, P is the gas pressure, p is the mass density, and oq is a 
positive constant. All of the physical quantities denote vertical averages at a given cylindrical radius 
r. The assumed variations of E, I, a, and 12 are the most natural extensions of the scalings used 
by BB into the pseudo-Newtonian regime. In § 3 we shall see that these scalings impose definite 
constraints on the mechanism responsible for producing the energy and mass loss that powers the 
outflow. 



2.2. Disk- Wind Coupling 

Our ultimate goal is to understand the possible role of the shear acceleration mechanism 
discussed by SBK in launching ADIOS-type outflows. Consequently, we focus here on gasdynamical 
outflows (winds) that exert no local torque on the disk, and therefore we stipulate that 12 = 12jy, 
where 

r 2 dM r 2 dM/dr 
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denotes the angular velocity of the wind at the point of contact with the disk. Substituting for I 
using equation (2.3) then yields 

which immediately implies that the variation of M is given by 

Ma(r 2 fi) AM , (2.7) 

or, using equation (2.2) for Q(r), 

A/(n -A) 



r 



3/2 



Moc _ . (2.8) 



We can relate the variation of M to the other physical quantities in the disk by writing 

M = 47rrHvp, (2.9) 

where H is the half-thickness of the disk and v is the radial velocity, defined to be positive for 
inflow. The value of H can be characterized by employing the equation of vertical hydrostatic 
equilibrium to write 

---7-^ W £ ' ( 2 - 10 ) 
z (r — 2) z r 

which is valid provided z <C r. A scale analysis of this expression yields the usual approximation 
(e.g., Narayan et al. 1997) 

which reduces to equation (20) of BB in the limit r — > oo. Since the disk must be thin for self- 
consistency, we require that ao < 1. 



2.3. Radial Momentum Conservation 

Our self-similar approach is based on the scalings for E, J, and a given by equations (2.3) and 
(2.4), respectively. We conjecture that the associated radial velocity v scales as 

v(r) = , (2.12) 

v< — 2 

where vq is a positive constant, implying that the radial velocity is proportional to the local freefall 
velocity, vg(r) = \/^/{r — 2). We shall verify this hypothesis by confirming that the resulting 
variation of the energy transport rate E(r) satisfies equation (2.3). Note that in order to avoid 



- 7- 



inflow faster than freefall, we require that Vq < 2. In a steady state, the radial equation of motion 
(eq. [6] of BB) can be written as 

dv ^ 9 1 1 dP 

v *- ra -irw-7* ■ (2 ' 13> 

We can obtain a differential equation for the pressure P by substituting for Q, p = P/a 2 , and v in 
equation (2.13) using equations (2.2), (2.4), and (2.12), respectively, which yields 



dlnP _ wg + 2 — 2 
dr ~ 2al(r - 2) 

with solution 



(2.14) 



P(r) = P (r - 2)K 2 +2^-2)/(2ag) ; (2 J5) 
where Po is a positive constant. The associated density variation is given by 

p(r) = P0 (r- 2 )K+ 2 ^+2^-2)/(2ag) ; (2 16) 

where 

PO = 4 ■ (2-17) 

Using equations (2.11), (2.12), and (2.16) to substitute for H, v, and p, respectively, in equa- 
tion (2.9) for the accretion rate, we obtain 

M = M r 3 / 2 (r - 2)("o+ 2 ng+2ag-2)/(2ag) ? ^.18) 

where 

M = 47ra « Po • (2.19) 

The power-law indices in this expression must match those appearing in equation (2.8), and there- 
fore we require that 

A = ^, vl = 2-2^l-Aal , (2.20) 

which together imply that the self-similar variation of the accretion rate is given by the unique 
solution 

r 3/2 

M{r) = M —— . (2.21) 

This solution for M(r) has an interesting behavior, in that it possesses a minimum at r = r ms = 6, 
which is the radius of marginal stability for a Schwarzschild black hole. The rate at which matter 
is fed into the wind therefore approaches zero as r — > r ms . This is physically reasonable since below 
this radius we expect that most of the remaining material will be swept directly into the black hole 
rather than expelled into the wind. For r > 6, M decreases monotonically with decreasing radius 
in response to the loss of mass into the wind, and as r — > oo, we find that M oc r 1 / 2 , in agreement 
with BB provided their parameter p = dlnM/dlnr = 1/2. Our sub-Keplerian assumption S7q < 1 
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can be combined with the requirement that Vq > to conclude that Og < (1 — Oq)/2, confirming 
that the disk is thin (see eq. [2.11]). Incorporating equations (2.20) into our previous results, we 
conclude that the pressure and density vary as 

p < r > - V^W • p(r) = • (2 - 22) 

which both increase monotonically with decreasing radius as required. In § 5 we discuss the possible 
relevance of this density distribution in accounting for the observations of time lags between the 
soft and hard X-rays in accreting galactic sources. 



2.4. Energy and Angular Momentum Conservation 

Following BB, we express conservation of angular momentum in the disk by writing the 
inwardly-directed angular momentum transport rate as 

J(r) = Mr 2 Q-g , (2.23) 

where Q{r) is the torque exerted by the disk interior to radius r on the outer portion of the disk. 
The outwardly-directed energy transport rate in the disk is likewise given by 

E{r)=GV-M ^-r 2 tf + \v 2 + h -a 2 -^-^j , (2.24) 

which incorporates the pseudo-Newtonian gravitational potential. Note that equation (2.24) in- 
cludes a term describing the radial kinetic energy of the gas, whereas BB ignored this term in their 
energy transport equation. The gas in the disk is assumed to have adiabatic index 7 = 5/3. We can 
solve for the torque by combining equations (2.3) and (2.23) and setting A = f^o/2 in compliance 
with equation (2.20), which yields 

Q(r) = I(r) = ^Mr 2 n . (2.25) 

The physical implications of this torque can be investigated by combining equation (2.25) with the 
expressions for M and f2 to determine the implied variation of the dynamic viscosity v. We discuss 
this issue further in § 3. By using equation (2.25) to eliminate g in equation (2.24) and substituting 
for 17 and a using equations (2.2) and (2.4), respectively, we can reexpress the energy transport 
rate as 

E(r) = e^±, (2.26) 

where 

e = l -\ a °-\ V « ■ (2 - 27) 
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Equations (2.20) and (2.27) establish three connections between the five parameters A, Qq, e, ao, 
and vq. Hence knowledge of any two is sufficient to calculate the remaining three. We shall therefore 
treat J7o and e as free parameters in our model and compute A, ao, and vq using the relations 

A = ^> a§ = 2ng-2e, vl = 8e - 10 n§ + 2 , (2.28) 

obtained using equations (2.20) and (2.27). 

The energy transport rate given by equation (2.26) does satisfy equation (2.3), and therefore 
we have confirmed the validity of our solution for v(r) given by equation (2.12). The establishment 
of the self-similar form for E stems from the fact that in our model, the local torque is equal to 
the local angular momentum transport rate (see eq. [2.25]), and consequently the outward energy 
transport due to the torque balances the inward energy transport due to the accretion of azimuthal 
kinetic energy. BB obtain the same result for gasdynamical winds when p = 1/2. The net energy 
transport is in the outward direction provided e > 0. As r — > oo, we find that E(r) — ► eMr~ l in 
agreement with the ADIOS model of BB. Substituting for M using equation (2.21), our results for 
the energy and angular momentum transport rates can be rewritten as 

3/2 3 3 

£ M=^o^, 7(r) = AM > 5(r)=AM„^. (2.29) 

Note that the torque Q possesses a minimum at the radius of marginal stability r = r ms = 6, 
but it does not vanish there. This behavior is consistent with simulations performed by Hawley 
& Krolik (2000), who find that the viscous stress is finite at r = r ms . Moreover, Narayan et al. 
(1997) suggest that the stress actually vanishes at the sonic point (between the horizon and the 
radius of marginal stability), where the radial velocity becomes supersonic and the gas therefore 
loses pressure support. It follows that the torque is finite at r = r ms , as implied by equation (2.29). 
However, it must be pointed out that this issue is a subject of ongoing debate. While Hawley & 
Krolik (2000), Agol & Krolik (2000), Gammie (1999) and Krolik (1999) contend that there is a 
finite stress at the radius of marginal stability, Armitage et al. (2000) and Paczynski (2000) argue 
that this should not be the case, at least for thin disks. In § 3 we examine the implications of 
our self-similar model for the variation of the viscosity and for the associated dissipation of kinetic 
energy in the disk. 



3. Viscosity and Dissipation 

Equation (2.25) for the torque Q in our self-similar RADIOS model has definite implications 
for the variation of the kinematic viscosity v. We can calculate the viscosity by writing the shear 
stress as (e.g., Frank, King, & Raine 1985) 

pur— , (3.1) 



^irr 2 H dr 



-10- 



and eliminating Q using equation (2.25), which yields 

" = -—{-ct) ■ (3 - 2) 

Substituting for and v using equations (2.2) and (2.12), respectively, we obtain 

_ v r(r -2) 1 / 2 



3r - 2 



(3.3) 



We can evaluate the variation of v most straightforwardly by calculating an equivalent value for 
the associated Shakura-Sunyaev a-parameter, defined by writing the shear stress as 

aP = - pur ^. (3.4) 
dr 

Setting P = pa 2 and combining the result with equations (2.2), (2.4), (2.28), and (3.3) yields 

( r \ 1/2 

a = a ( — - j , (3.5) 



where 



_ fi (8 e- 10 nl + 2) 



1/2 



a " s 4 ( ng-°) • <3 ' 6) 

We plot cko as a function of e and f2o m Figure 1. Equation (3.5) indicates that a has a very weak 
radial dependence, increasing by a factor of 1.2 as r decreases from infinity to r = r ms = 6. Hence 
the a-parameter associated with our self-similar RADIOS model is practically constant within the 
hot inner region of the disk, which is consistent with the conventional assumption that the shear 
stress is approximately proportional to the pressure (e.g., Narayan et al. 1997). Based upon the 
observed behavior of a, we conclude that our model has an implied viscosity variation that is 
consistent with fundamental notions regarding the transport of angular momentum in the shear 
flow. This tends to further support our fundamental assumptions regarding the scalings of E, I, 
and a. 



3.1. Implications for Dissipation and Cooling 

Next we consider the implications of our self-similar model for the dissipation on energy in the 
shear flow, and for the rate of cooling from the disk due to the loss of energetic particles into the 
wind. The equation governing the internal energy of the matter in the disk is 

jj^ — 1 jj^ ' U viscous L'wind j ("•') 

where D/Dt represents the comoving (Lagrangian) time derivative following the gas, 7 is the 
adiabatic index, U = P/("y — 1) is the internal energy density, and f/ v i SC ous and U w - m d denote the 
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(positive) rates of viscous heating and cooling into the wind, respectively. Using the standard 
expression for the viscous heating rate (Frank, King, & Raine 1985), 



(Adscous = pvr [ — 1 , (3.8) 



we obtain in a steady state 



vp da 2 2 dp 2 / dn\ 

+ va — = pur — -U wind . (3.9) 



1 — 7 dr dr \dr 

This is equivalent to equation (2.12) from Narayan et al. (1997), except that we have chosen v 
to be positive for inflow. In the standard ADAF theory, a factor / is introduced to represent the 
fraction of the energy dissipated via viscosity that goes into heating the gas, with the remainder 
escaping as radiative losses. We can define an analogous heating efficiency factor / in the context 
of our disk/ wind model by writing 

( dQ\ 2 2 ( dQ s 



s » vr =""• I* J - -" (3 ' 10> 

Hence 1 — / gives the fraction of the energy dissipated via viscosity that is expelled into the wind 
via the internal energy of the escaping particles. Narayan et al. and other authors focusing on the 
standard ADAF scenario usually assume that / = 1 throughout the disk. However, in our RADIOS 
model we can calculate / self-consistently by solving for f/ w i n d and then utilizing equation (3.10). 
Setting 7 = 5/3 and substituting for v, a 2 , p, v, and Q using our self-similar model, we find that 
equation (3.9) yields for the implied wind cooling rate 

"fig (6-r) + 4e(r-2)" 



j r _ Po vq 



(3.11) 



(r - 2) 9 / 2 

The corresponding result for the heating efficiency factor obtained using equation (3.10) is 

™- 

Since the energy dissipated from the shear flow via viscosity is divided between heating of the 
thermal particles and acceleration of the nonthermal (wind) particles, it follows that / must lie in 
the range < / < 1 if we are to obtain a physically reasonable solution. In order to satisfy this 
condition as r — > oo, we must stipulate that 

e > -f ■ (3-13) 

Note that the heating efficiency / displays a weak radial dependence, increasing by a factor of 4/3 
as r varies from r = 6 to infinity. This indicates that a nearly fixed fraction of the energy supplied 
via the dissipation of the shear goes into heating of the disk, with the remainder powering the 
outflow. This seems to be a reasonable result, and one cannot say more in the absence of a detailed 
model for the acceleration mechanism that energizes the wind particles. We plan to address this 
issue by developing an advection-dominated disk/wind model that incorporates the Fermi-shear 
acceleration process discussed by SBK in future work. 
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3.2. Parameter Constraints 

We have assumed that the disk is thin and sub-Keplerian, and that the inflow velocity is less 
than the local freefall velocity. These assumptions can be combined to constrain the range of values 
for the various parameters appearing in the model. In order for the disk to be thin, we require 
that <2q < 1. Equation (2.28) for Oq then implies that for a given value of Clo, the parameter e must 
satisfy 

Max (o , nl-^j < e < n 2 , (3.14) 

where the upper bound follows from the requirement that Oq > 0. By combining equation (3.14) 
with our sub-Keplerian restriction O < 1 , we find that e cannot exceed unity in general. We must 
also stipulate that < Vq < 2 in order to keep v 2 positive while avoiding inflow faster than the 
free-fall velocity. Equation (2.28) for Vq then implies that e must fall within the range 

Max f , 1 J < e < . (3.15) 

We plot the constraints represented by equations (3.14) and (3.15) in the (CIq , e) parameter space 
in Figure 2. The intersection of the constraints is expressed by the condition 

Max (o , 5 J3_Z± \ <e (3.16) 

We can combine this condition with equation (3.13) to conclude that in order for our model to be 
fully self-consistent, e and must satisfy 

Max (^L , <e<Q 2 , < fig < 1 . (3.17) 

The allowed region of the (S7o , e) parameter space is the curved, shaded region depicted in Figure 2. 
For convenience, we have also included the shaded region of self-consistency in Figure 1 for a®. Note 
that we can obtain values for «o in the range < «o < 1 within the region of the parameter space 
satisfying the conditions given by equations (3.17). In § 4 we present a quantitative discussion of 
the implications of the pseudo-Newtonian potential for the coupled disk/wind structure. 



4. EFFECTS OF THE RELATIVISTIC POTENTIAL 

In this section we compare the results obtained using our self-similar, pseudo-Newtonian, 
relativistic advection-dominated inflow-outflow solution (RADIOS) with the ADIOS results derived 
by BB using the standard Newtonian potential. Before proceeding, it is important to point out 
that the ADIOS model includes the possibility of a mismatch between £1 and £lwi i- e -> there may 
be a torque between the disk and the wind. Since the winds treated here are gasdynamical, we 



-13- 



have assumed that no such torque exists, and consequently we have set f2 = This corresponds 
to case (iv) of BB, with A = A cr i t , where 

(lOe + Ap - Aep) 1 1/2 



A crit = 2p 



(4.1) 



_(2p+l)(4p 2 + 8p + 15)_ 

Furthermore, as we have already pointed out, our utilization of the pseudo-Newtonian potential 
imposes a length scale (i.e., the horizon at r = 2) that implies a definite behavior for M as given by 
equation (2.21). This in turn requires that we set p = 1/2 in the ADIOS model in order to match 
the behavior of our M solution as r — ► oo. Incorporating this value for p into equation (4.1) then 
yields 

1 /l +4e \V2 

A =2-(^J ' ^ 
which can be combined with equations (19) and (20) of BB to obtain 

/l + 4e\ 1/2 /2-2e\ 1/2 , A N 

These results are consistent with our equations (2.28) in the special case vq = 0, which is due to the 
fact that BB neglected the radial kinetic energy term in their energy transport equation. We shall 
use the subscript "^4" to denote quantities related to the ADIOS model. Unsubscripted quantities 
will continue to refer to our RADIOS model. The self-similar ADIOS expressions for the angular 
velocity and the isothermal spound speed are then given by 

Q A {r) = n r" 3 / 2 , a A {r) = a r' 1 ' 2 , (4.4) 

respectively. The corresponding ADIOS solutions for the energy, angular momentum, and mass 
transport rates reduce when p = 1/2 to 

E A (r) EeMpr 1 / 2 , I A (r) = XM r, M A (r) = M r 1 ' 2 , (4.5) 

respectively. These solutions obviously agree with our equations (2.21) and (2.29) in the limit 
r — > oo as required. In Figures 3 and 4 we compare our results for the energy and mass transport 
rates E{r) and M(r) with the corresponding ADIOS quantities E A (r) and M A (r) for radii in the 
range r ms < r < 10, where r ms = 6. Note that for the same value of e, energy is transported 
outward through the disk at a substantially greater rate in the relativistic (pseudo-Newtonian) 
model due to the steepness of the potential in the inner region of the disk. 



4.1. Disk/Wind Structure 

The modifications to the ADIOS disk structure due to the incorporation of the pseudo- 
Newtonian gravitational potential in RADIOS can be summarized by writing 

M_ f^\ 2 (I_V' 2 (E_\ 1/2 



M A n A \a A J \I A J \E A ) r-2' (4 ' 6) 
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which together suggest that no dramatic changes in the disk properties occur until one approaches 
the horizon at r = 2. In fact, the importance of the subtle variations in structure between the 
relativistic and nonrelativistic models really becomes clear only when one examines the properties 
of the associated outflows. Due to the symmetry across the disk midplane, identical outflows 
emanate from the upper and lower surfaces of the disk. For brevity, we shall refer to these two 
outflows collectively as the "wind." The total amounts of mass and energy emitted into the wind 
per unit time between radii r ms and r are given by 

m(r) = M(r) - M(r ms ) , L(r) = E(r ms ) - E(r) , (4.7) 

for the pseudo-Newtonian model. The corresponding expressions in the Newtonian case are 

m A {r) = M A (r) - M A (r ms ) , L A {r) = E A (r ms ) - E A (r) . (4.8) 

In the absence of any outflows (e = 0), the energy and momentum transport rates are conserved 
and therefore L(r) = L A (r) = 0, as in the ADAF models of Narayan & Yi (1994). At large radii, 
we find that the total power emitted into the wind approaches 

• 6 3 / 2 

L(oo) = E(r ms ) = eM — , (4.9) 

and 

L A (oo) = E A (r ms ) = eM 6- 1/2 , (4.10) 

for the pseudo-Newtonian and Newtonian models, respectively. If we interpret the radius of 
marginal stability r ms = 6 as the inner edge of the disk, then the values L(oo) and L^oo) represent 
the total power emerging from the disk into the wind from all radii r > 6. Equations (4.9) and 
(4.10) indicate that L(oo)/L A (oo) = 2.25, so that the relativistic potential roughly doubles the total 
power in the outflow for fixed e and Mq. In an actual accretion disk, we expect that the wind will 
terminate at some outer radius r out > r at which the hot inner region transitions into the optically 
thick, cool outer disk. We are not particularly concerned with the value of r out here, although we 
note that most hot disk models give values in the range r out ~ 10 2-4 (Narayan, McClintock, & Yi 
1996; Narayan, Barret, k McClintock 1997). 

It is also interesting to compute the amount of mass and energy emitted into the wind per unit 
radius per unit time. Utilizing our pseudo-Newtonian RADIOS expressions yields 

dL ■ r x / 2 (r + 6) dm ■ r 1 / 2 (r - 6) , A „. 

The corresponding expressions obtained using the standard ADIOS model of BB are 

dL A • r~ 3 / 2 dm A • r -1 / 2 

*- = eM °— ■ ir-*>— ■ (412) 

These results are compared graphically in Figures 5 and 6 for r in the range 6 < r < 10. For a 
given value of e, the power emitted into the wind is significantly greater in the pseudo-Newtonian 
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case than in the standard Newtonian model due to the steeper potential gradient. Note that in 
the pseudo-Newtonian case, dm/dr = at the radius of marginal stability r ms = 6. As explained 
earlier, we therefore interpret r ms as the inner edge of the disk, and assume that essentially all of 
the matter remaining at that radius falls directly into the black hole. 



4.2. Bernoulli Parameters 

Another useful way of contrasting the properties of the Newtonian and pseudo-Newtonian 
inflow-outflow models is to calculate the specific energy (i.e., the Bernoulli parameter) in the disk 
and the wind as functions of radius. In the pseudo-Newtonian case, the disk Bernoulli parameter 
is given by 

B{r) = -v- + -r + - a- - — = ^-^ °2- , (4.13) 

where we have utilized equations (2.28) to arrive at the final result. The corresponding expression 
in the Newtonian case is 

obtained using equations (4.3) and (4.4). In equations (4.13) and (4.14), r is interpreted as the 
cylindrical "starting radius" at the base of the wind in the disk. After leaving the disk, the wind 
is expected to expand to larger cylindrical radii as the outflow proceeds. Note that equation (4.13) 
for B includes the contribution due to the radial kinetic energy, whereas the expression for B& 
employed by BB (eq. [4.14]) assumes that v$ = 0. 

In general, the pseudo-Newtonian disk Bernoulli parameter B may be either positive or negative 
depending on the values of e, Q,q, and r. For given values of Qq and r, we find that B < if and 
only if e < e*, where 

_ 4-2r + (9r-20)ng 
£ * = fo~l2 ' (4 ' 15j 

The roots of B occur at radius r = , where 

r. = 12£ - 2 °^ + 2 4 , (4.16) 

with B < for r > r* and B > for r < . The situation is complicated by the fact that e and f^o 
must also satisfy the self-consistency conditions given by equations (3.17). However, we can make 
a few general statements; (i) If $7q < 1/3, then B > for all radii r > 6; (ii) If 1/3 < Sl§ < 4/5, 
then there is an allowed range of e values for which B changes sign at r = r*, where r* > 6; (iii) If 
> 4/5, then B < for all radii r > 6. Hence, in the pseudo-Newtonian case, the inner region of 
the disk may be unbound even when the outer disk is bound. 

In Figure 7 we plot B and Ba as functions of radius for various values of e, with vo = in 
order to facilitate a direct comparison between the RADIOS and ADIOS model results. Recall 
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that setting vq = implies that f2o = \/(4e + l)/5 according to the last of equations (2.28). In 
the Newtonian case, we find that Ba > for all values of r if e < 1/6, and Ba < for all radii if 
e > 1/6. In the pseudo-Newtonian case with vo = 0, we also find that B > for all values of r if 
e < 1/6, but if e > 1/6, then B changes sign at = 20e/(6e — 1). Note that B exceeds Ba at all 
radii, which is due to the effect of the pseudo-Newtonian potential. 

A positive disk Bernoulli parameter violates self-consistency in standard ADAF models because 
these do not include outflows. By contrast, in the ADIOS and RADIOS models considered here, 
outflows are included and therefore we can obtain self-consistent solutions even when the Bernoulli 
parameter is positive, although the most plausible scenario is one in which it is driven to a negative 
value as a result of the outflow. Specifically, the most energetic outflows are obtained when e = 1, 
which yields a negative (bound) value for the disk Bernoulli parameter for all values of r. Conversely, 
the most positive values for the Bernoulli parameter are obtained when e = 0, which corresponds 
to the no-wind solution since in this case the energy transport rate in the disk vanishes. The source 
of the change in behavior of the RADIOS results plotted in Figure 7 relative to those obtained 
using the ADIOS model can be traced to the quasi-Keplerian form for O operative in the pseudo- 
Newtonian case, expressed by equation (2.2). As a consequence of equation (2.2), the azimuthal 
kinetic energy term (l/2)r 2 fl 2 in equation (4.13) for B varies as r/(r — 2) 2 , whereas the terms 
describing the internal and gravitational potential energies each vary as l/(r — 2). At sufficiently 
small radii, the azimuthal kinetic energy (which is positive) therefore overwhelms the remaining 
terms, which have the weaker l/(r — 2) radial dependence. The vigorous azimuthal motion is due 
to the stronger gradient in the pseudo-Newtonian potential relative to the Newtonian potential. In 
the ADIOS model all three terms in equation (4.14) for Ba have the same scaling, and therefore 
Ba has the same sign at all radii. 

The Bernoulli parameter for the wind particles is given in the RADIOS model by 

dL = _dE / d L= «(r + 6) 
V ' dm dM/dr (r - 6) (r - 2) v ; 

The equivalent expression obtained using the ADIOS model results is 

Mr) - P- = ~ dEA,dr = ~ ■ (4.18) 
dm A dM A /dr r 

Note that b and bA are always positive, and that they each exceed the associated disk Bernoulli 
parameters B and Ba, respectively. This implies that the energy per unit mass in the wind is 
higher than that in the disk, which is consistent with our disk/wind scenario. As the wind flows 
away from the black hole, we expect the gas to expand adiabatically, passing through a critical 
point at which the flow velocity equals the sound speed as discussed by SBK. In this situation, the 
specific energy of the wind particles is conserved and resides ultimately in the kinetic energy of 
the outflow. It must be emphasized that the wind Bernoulli parameter describes only the kinetic, 
potential, and internal energy per unit mass, since it is based on the energy transport equation in 
the disk, which is classical in form (see eq. [2.24]). Hence the rest mass energy of the wind particles 
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is not included in either b or b A . It follows that if the outflowing gas expands adiabatically, then 
the total energy per unit mass in the cold, terminal wind is given in the pseudo-Newtonian and 
Newtonian cases by (e.g., SBK) 

r«6 + l, T A ^b A + l, (4.19) 

respectively, where T and F A are the corresponding asymptotic Lorentz factors of the outflow in 
the two models. 

Without considering the structure of the outflow in any detail, we can use equation (4.19) to 
determine the asymptotic nature of the wind originating at radius r. In the pseudo-Newtonian 
RADIOS case we can combine equations (4.17) and (4.19) to obtain 

r(r )a ,Kr) + l= (r ^+ g > 2) + l. (4.20) 

In the Newtonian ADIOS case the corresponding result obtained using equations (4.18) and (4.19) 
is 

T A (r) « b A (r) + 1 = 1 + 1. (4.21) 

According to equations (3.17), the maximum value of e satisfying our self-consistency conditions 
is e = 1, corresponding to the limiting case f2o = 1- We find that in the pseudo-Newtonian 
RADIOS model with e = 1, the asymptotic Lorentz factor r > 2 for r < 8.3, and r > 3 for r < 7.3. 
This indicates that the inner region of the disk can produce a relativistic wind/outflow that may 
subsequently transition into a relativistic jet at large distances from the black hole, provided the 
gas expands adiabatically. Conversely, in the corresponding ADIOS model with e = 1, the largest 
possible value for T A is 1.17 (obtained for starting radius r = r ms = 6), and therefore the standard 
ADIOS outflow is clearly nonrelativistic. These results are depicted in Figure 8. 



5. DISCUSSION AND CONCLUSION 

Our relativistic, advection-dominated inflow-outflow solution (RADIOS) incorporates a pseudo- 
Newtonian gravitational potential, and therefore it provides an approximate representation of the 
effects of general relativity on the properties of the accretion disk and the associated outflow, which 
is assumed here to be a gasdynamical wind. Consideration of mass, momentum, and energy con- 
servation yields a new, self-similar solution for the disk/wind structure that is rather different from 
that obtained by BB. In particular, our results show that a powerful relativistic outflow originates 
from the inner region of the disk, close to the radius of marginal stability. The appearance of this 
outflow is a consequence of the steep gradient in the pseudo-Newtonian potential. Consequently, the 
RADIOS model may provide a better description of the inner region of the disk than the standard 
ADIOS model, which is based on Newtonian gravitation. Although it utilizes a pseudo-Newtonian 
potential, our model nonetheless incorporates the same fundamental physical scaling principles as 
the ADIOS model. Namely, we assume that the angular velocity Q scales as the Keplerian value 
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Uk (see eqs. [2.1] and [2.2]), and we assume that the energy transport rate per unit mass E/M and 
the specific internal energy a 2 each scale as the gravitational potential energy <I> = —l/(r — 2). We 
also assume that the angular momentum transport rate per unit mass I/M scales as the specific 
angular momentum r 2 fl. 

Evidence supporting our assumption regarding the variation of f2 is provided by the simulations 
of Narayan et al. (1997), who recover this behavior over several orders of magnitude in radius 
(outside the radius of marginal stability) for a ~ 0.1 — 0.3. Although Narayan et al. (1997) assume 
that there is no outflow from the disk, we can still utilize our equation (2.2) based on their work 
because we assume here that there is no relative torque between the disk and the wind. Hence we 
expect little change in the angular momentum distribution as a result of the outflow. This also 
supports our assumption regarding the variation of the angular momentum transport rate / since 
in a self-similar model such as ours, the scaling of / follows directly from the variation of f2 by 
virtue of equation (2.23). 

The pseudo-Newtonian potential introduces a length scale into the problem due to the existence 
of the horizon at r = 2. No such length scale appears in the Newtonian ADIOS model of BB. The 
imposition of a length scale in our model results in a definite behavior for the variation of the 
accretion rate given by equation (2.21), which approaches M oc r 1 / 2 as r — > oo. When coupled 
with our assumption of a torque-free (gasdynamical) wind, the number of free parameters in our 
model is reduced to two, which we take to be f^o an d e. The remaining constants A, cto, and vq 
can be computed in terms of f^o an d e using equations (2.28). The self-similar solutions we obtain 
asymptotically approach the corresponding, Newtonian ADIOS results for a gasdynamical wind as 
r — > oo if we set p = din M/ din r = 1/2 and A = A cr i t in the ADIOS model (see eq. [4.1]). Hence the 
gravitational physics near the horizon forces the flow structure to satisfy M oc r 1 / 2 at large radii. 
In the Newtonian ADIOS model satisfying this constraint, the outflow is always nonrelativistic. In 
our model, the outflow can be highly relativistic, implying a possible connection with large-scale 
jets. 

The asymptotic Lorentz factor in the wind, F, formally diverges as r — > 6 because dm/dr 
vanishes at the radius of marginal stability r ms = 6. This is due to the fact that for r < 6, the disk 
is unstable and all of the remaining material is swept directly into the black hole (see, however, 
our discussion of the implications of transonic flow below). Conversely, the differential power in 
the outflow, dL/dr, remains finite even at r = 6 because energy continues to be transported in 
the outward radial direction by the self-similar torque, which has a radial dependence given by 
equation (2.25). Even with our self-consistency conditions (eqs. [3.17]) taken into consideration, 
equation (4.20) indicates that the asymptotic Lorentz factor for the outflow F can still greatly exceed 
unity for the portion of the wind emitted close to the radius of marginal stability. By contrast, 
when the same parameters are used in the ADIOS model of BB, the outflow is subrelativistic from 
all radii in the disk (cf. eq. [4.21]). 

The radial Mach number is conserved in our self-similar model, and therefore the flow remains 
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subsonic in the entire computational domain (r > r ms ). Of course, in reality the gas must certainly 
achieve a supersonic velocity before it crosses the event horizon. The transition to supersonic flow 
requires that the material pass smoothly through a sonic point, which is a critical point for the flow. 
The existence of the sonic point therefore imposes critical conditions that the overall dynamical 
structure must satisfy. This does not appear to be a major problem for our model because the sonic 
point is expected to lie between r = r ms and the horizon r = 2, which is outside the computational 
domain of our disk model. However, the issue can be addressed quantitatively by constructing a 
global model in which the subsonic region (r > r ms ) is described by our self-similar RADIOS solution 
and the transonic region (2 < r < r ms ) is described by the solution of Narayan et al. (1997; see 
also Chen et al. 1997). A more sophisticated alternative would be to modify the Narayan et al. 
calculation to include mass loss such as that envisioned here. We expect that in the resulting 
model, matter will actually be expelled into the wind all the way down to the sonic radius (located 
between the radius of marginal stability and the horizon), where the pressure support for the wind 
essentially vanishes. Hence, there should still be some minimal amount of outflow at r = r ms , and 
consequently T should remain finite at the radius of marginal stability, in contrast to our results. 
Nonetheless, we expect that the dynamical structure of such a generalized model will be similar 
to that of the RADIOS model because the Q(r) distribution that we utilize is consistent with that 
obtained by Narayan et al. We therefore believe that our model provides an intriguing glimpse 
into the possible behavior to be expected when fully relativistic calculations that include particle 
acceleration and mass loss become available. 



5.1. Shear Acceleration and the Origin of the Outflows 

In the treaments of advection dominated inflow-outflow solutions that have appeared to date 
in the literature (BB; Donea et al. 1999), the existence of outflows is generally assumed without 
any elaboration on the physical mechanisms that could give rise to them. That is also the case in 
the work presented here, since the outflow power is calculated by differentiating the disk energy 
transport rate E, without making any reference to a microphysical model. Although several mech- 
anisms for producing outflows have been proposed (see references in SBK), few make a connection 
between the outflow and the underlying accretion disk. Das & Chakrabarti (1999) and Das (1999) 
describe pressure-driven winds from centrifugally-supported boundary layers and shocks in the in- 
ner regions of disks. Xu & Chen (1997) proposed an advection-dominated flow where the central 
object effectively acts as a scatterer which redirects the inward flow at low latitudes into an outflow 
at high latitudes. The above-mentioned models compute the mass outflow rates self-consistently, 
but the associated winds are subrelativistic and therefore the possible connection with large-scale 
relativistic jets is unclear. 

The formation of a relativistic outflow requires the operation of a mechanism in the disk that 
preferentially accelerates high energy particles. The resulting particle distribution comprises a 
thermal component that is bound to the disk and a nonthermal component that is unbound, repre- 
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senting the particles that are destined to escape into the wind. SBK have described a second-order 
Fermi acceleration mechanism powered by the shearing magnetic field that is capable of producing 
a relativistic outflow. Some related conceptual development has also been carried out by Blackman 
(1999) and by Gruzinov & Quataert (1999). The physical scenario in which the acceleration of 
the wind protons takes place is the same as that in which hybrid viscosity (Subramanian, Becker, 
& Kafatos 1996) is operative. The azimuthal kinetic energy in the shear flow is dissipated partly 
by viscous heating of the thermal protons, and partly by second-order Fermi acceleration of the 
nonthermal protons, which in turn form a relativistic outflow. SBK demonstrated that this mech- 
anism can lead to the energization of protons via shear acceleration in a tenuous corona overlying 
a hot accretion disk. The resulting specific energy in the corona can be large enough to drive a 
relativistic outflow. In their work, SBK focused on shear acceleration occurring in a tenuous corona 
because the density in the underlying (luminous) disk is so large that collisional losses overwhelm 
shear acceleration there. However, in the advection-dominated accretion flows envisioned here, the 
low density environment of the (underluminous) disk itself is expected to be a favorable site for 
this acceleration mechanism. 



5.2. Model Predictions 

The unique solution for the M variation given by our equation (2.21) and the associated density 
distribution p oc (r — 2) _1 have definite consequences for the timing properties of accretion-powered 
sources. It has been proposed (Kazanas, Hua & Titarchuk 1997; Hua, Kazanas & Titarchuk 1997; 
Hua, Kazanas, & Cui 1999; Kazanas & Hua 1999) that the timing properties of these sources 
(light curves, power spectra, coherence, and time lags between soft and hard X-rays) can be easily 
understood as the result of Comptonization of soft photons produced in the vicinity of the com- 
pact object by a hot (10 s K) corona that extends to r ~ 10 3 Schwarzschild radii. In particular, 
it was pointed out that such an interpretation provides a means for probing the density structure 
of the Comptonizing corona, mainly as a consequence of the dependence of the X-ray hard lags 
on the Fourier frequency. Hua, Kazanas, & Cui (1999) analyzed and fit the lags from a number 
of accreting galactic black hole candidate sources. In a few cases, the lag data were consistent 
with the density distribution p oc r -3 / 2 implied by the ADAF model (Narayan & Yi 1994). How- 
ever, in most cases the lag data were consistent with a corona density profile p(r) oc r _1 , which 
is the behavior that is uniquely determined by our self-similar RADIOS model for r>2. Hence 
our pseudo-Newtonian, self-similar RADIOS model may provide the first explanation for a ubiqui- 
tous observational conclusion regarding the apparent radial density variation in accreting compact 
sources. 
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5.3. Conclusion 

We have described the overall framework of a pseudo-Newtonian RADIOS model that asso- 
ciates relativistic outflows with advection-dominated accretion disks. The self-similar solutions 
obtained for the physical properties of the disk extend over several decades in radius and therefore 
provide a meaningful description of the global disk structure. By extending our model to large 
radii, we are able to constrain the parameters in the associated Newtonian ADIOS model. Our 
analysis yields unique radial dependences for the accretion rate, the density, and the pressure. For 
the special case of a torque-free, gasdynamical wind, which is the focus here, we find that the 
density varies as p oc r _1 as r — ► oo. This density variation is consistent with that implied by 
observations of X-ray time lags in accreting galactic sources. Moreover, the Shakura-Sunyaev a- 
parameter implied by our model remains nearly constant, independent of the viscosity model (see 
cq. [3.5]). This is consistent with the basic notion that the shear stress should be roughly propor- 
tional to the pressure. When combined with the fact that our assumed variation of Q is consistent 
with published numerical simulations, we conclude that our fundamental scaling assumptions are 
reasonable. Nonetheless, these assumptions cannot be fully justified a priori, and need to be tested 
by global calculations. In our RADIOS model, the strength of the outflow and the value of the 
mass loss index p = d In M/d In r are determined by the shape of the non-Newtonian potential close 
to the black hole. However, it is possible that the behavior p — > 1/2 as r — ► oo implied by observa- 
tions is determined, not by the shape of the potential, but rather by the details of the mechanism 
driving the outflow (e.g., second order Fermi acceleration, magneto-centrifugal forces, convection). 
If the latter is true, then the self-similar solution presented in this paper, though mathematically 
attractive, may have limited physical relevance. More work on this issue is needed. 

The RADIOS model we have developed can produce highly relativistic outflows, whereas the 
outflow is always nonrelativistic in the corresponding ADIOS model. Furthermore, like ADIOS, 
our model does not suffer from the strong sensitivity to the value of the adiabatic index 7 that is 
displayed by the standard ADAF models, which require 7 < 5/3 (e.g., Narayan & Yi 1994, 1995a). 
In contrast to the results for bipolar ADAF outflows obtained by Igumenshchev & Abramowicz 
(1999), which require a ~ 1, Figures 1 and 2 suggest that in our model outflows can occur for 
values of a well below unity. This work may therefore help to establish a physical connection 
between advection-dominated accretion flows and the relativistic outflows possibly emanating from 
underluminous black hole systems such as Sgr A* (e.g., Falcke 1999; Donea et al. 1999). The 
model may also be applicable to the mildly relativistic outflows observed from galactic black hole 
candidates (Mirabel & Rodriguez 1999; Rodriguez & Mirabel 1999). However, the nature of the 
fundamental acceleration mechanism that energizes the wind particles is still an open question. One 
possibility discussed here is second-order Fermi acceleration driven by the shear of the magnetic 
field lines. Additional work needs to be done in this area before the complete disk/wind/jet scenario 
can be adequately understood. 

The authors thank the anonymous referee for a variety of useful comments. PAB also gratefully 
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FIGURE CAPTIONS 

Fig. 1. - Contour plot of the Shakura-Sunyaev viscosity normalization parameter qq (eq. [3.6]) 
as a function of e and f^o- The contour values are ao = 1 (solid line), ao = 0.5 (dashed line), and 
ao = 0.01 (dotted line). In the region between the solid lines, ao exceeds unity. Also shown for 
reference is the shaded region of self-consistency for our model, given by equations (3.17). Note 
that this region contains values for ao in the physically acceptable range < ao < 1. 

Fig. 2. - Sections of the (Qq , e) parameter space corresponding to the various constraints given 
by eqs. (3.17). The region above the solid line satisfies e > J7q/4 , the region above the dashed line 
satisfies e > (5 J7q — 1)/4 , and the region below the dotted line satisfies e < Qq . The curved, shaded 
area is the intersection of these constraints and therefore it represents the region of self-consistency 
for our model, which is also indicated in Fig. 1. See the discussion in the text. 

Fig. 3. - Comparison of the pseudo-Newtonian and Newtonian energy transport rates E 
(eq. [2.29]; solid line) and Ea (eq. [4.5]; dashed line) in units of eMo as functions of the radius 
in the disk r. In the inner region, E exceeds Ea due to the steepness of the pseudo-Newtonian 
potential. The two results converge as r — ► oo. 

Fig. 4. - Pseudo-Newtonian and Newtonian accretion rates M (eq. [2.21]; solid line) and Ma 
(eq. [4.5]; dashed line) in units of Mq plotted as functions of the radius in the disk r. The accretion 
rate M has a minimum value at r = 6, which is the radius of marginal stability, and the inner edge 
of our disk model. At large radii, M approaches the standard ADIOS result given by Ma- 

Fig. 5. - Differential luminosities dL/dr (eq. [4.11]; solid line) and dLA/dr (eq. [4.12]; dashed 
line) emiited into the wind are plotted in units of eMo as functions of the starting radius in the disk 
r. These two functions correspond to the pseudo-Newtonian and Newtonian models, respectively. 
The power emitted into the wind is greatest is the pseudo-Newtonian case. 

Fig. 6. - Comparison of the pseudo-Newtonian and Newtonian differential mass-loss rates 
dm/dr (eq. [4.11]; solid line) and driiA/dr (eq. [4.12]; dashed line) emitted into the wind, plotted 
in units of Mq as functions of the starting radius in the disk r. In the pseudo-Newtonian case, this 
rate approaches zero close to the radius of marginal stability, and any remaining material is swept 
directly into the black hole. 

Fig. 7. - Pseudo-Newtonian and Newtonian disk Bernoulli parameters B(r) (eq. [4.13]; solid 
line) and Ba(t) (eq. [4.14]; dashed line) plotted as functions of the radius in the disk r for vq = 0. 
In panel (a), e = 0.1, and both B and Ba remain positive at all radii. In panel (b), e = 0.3, and 
B changes sign at r = 7.5, whereas Ba remains negative at all radii. This difference in behavior 
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is due to the strong variation of the quasi-Kcplerian angular velocity 0, in the inner region of the 
disk in the pseudo-Newtonian model. 

Fig. 8. - Asymptotic Lorentz factors T (eq. [4.20]; solid line) and Ta (eq. [4.21]; dashed line) 
for the outflows in the pseudo-Newtonian and Newtonian models, respectively, plotted as functions 
of the starting radius in the disk r. The Newtonian result is nonrelativistic for all values of r, 
whereas the pseudo-Newtonian result is strongly relativistic in the inner region of the disk, and 
formally diverges as r approaches the radius of marginal stability r ms = 6. 
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